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1 Introduction 


At low energies New Physics (NP) was expected to be seen indirectly in the down 
qnark sector. The LHC offered chance for direct search of NP. Althongh, NP is not 
found yet directly at high energies, there are a number of reasons why we still ex¬ 
pect to find its presence. For example on experimental side, in B physics tensions 
between SM expectations and experimental results are found. It was noticed that in 
B —)■ K*ij,~^ij,~ observable, known as P^, deviates for about 3a [T] from the SM pre¬ 
diction, the ratio = BR{B —)■ DrUr)/BR{B —)■ exhibit 3.8cr deviations 

[2] and Rk = BR{B —)■ K/BR{B —)■ A'e+e“)g 2 g[i^e]Gev 2 has 2.6cr dis¬ 
crepancy from the SM value |3]. Two of these observables in H iPh) transitions 
are result of flavor changing neutral current processes (FCNC), while the one for 
B L)h) is a result of charge current. These anomalous results stimulated numerous 
studies of NP in B meson system. One should keep in mind that all other B meson 
physics observables offer additional constraints to new physics. 

Top quark physics seems to be important for NP searches in the up-like quark sector. 
Properties and dynamics of top quarks attract a lot of attention on experimental and 
theoretical side. However, there is no indication yet about NP presence. The question 
one can ask: is there any chance to observe NP effects in charm FCNC physics? The 
constraints on NP in semileptonic charm decays driven by charge currents have been 
discussed in Ref. |1]. On the other hand, FCNC rare charm processes are accessible 
in radiative or semileptonic decays in which transitions c —>■ uy and c —)■ uPl~ occur. 
The main obstacle to search for NP in rare charm decay is the presence of many 
non-charm resonances in the vicinity of D mesons masses. Strong role of GIM mecha¬ 
nism is very important in charm FCNC dynamics. The interplay of CKM parameters 
and masses of down-like quarks leads to strong suppression of all FCNC in D meson 
processes. In addition, long distance contributions overshadow short distance effects. 
The main issue is how to separate information on short distance dynamics, either 
within SM or in its extensions. This is a longstanding problem in rare charm decays. 
Three years ago flavor community was concerned about discrepancy between mea¬ 
sured and expected CP violating asymmetry in charm decays [S]. Although this 
discrepancy seems to disappear, many studies and additional checks of the observed 
anomaly in rare charm decays were performed. The question on observability of CP 
violation in charm rare decays, should be answered. In Sec. 2 contributions to c —)■ uy 
and c — uPl~ decay modes are reviewed. The exclusive weak radiative D —)■ Hy 
decays are discussed in Sec. 3. and D —D —P(P')p+/i“ were analysed in 
Sec. 4. Tests of CP violation in charm meson decays with the leptons in the hnal 
state are discussed in Sec. 5. Last section contains the summary. 
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2 Inclusive decay modes: c ^ and c —)► ul^l 


The c ^ u'y and c —)■ ul^l transitions within SM can be approached by the effective 
low-energy Lagrangian: 


rSD _ 
^eff - 


4Gf 


V:tVub E 

i=7,9,10 


The operators are then: 


Qi = 

e 

■^mcFpyUapAl + 75) c , 

Q9 = 

6^ 

Qio = 



( 1 ) 


( 2 ) 


In Q Ci denote as usnal Wilson coefficients (they are determined at the scale /i = 
iTT'c), is the electromagnetic field strenght and ql = ^(1 ~ the case of 

c ^ U'j decay only Cy contribntes, while in the case of c —)■ ul^l~ all three Wilson 
coefficents are present. The QCD corrections enhance the rate to BR{c —>■ w))sm = 
2.5 X 10“® [SI E]. Within Standard model the short distance contribution coming 
from Qy ,9 leads to the branching ratio BR{D —)■ ~ 3.7 x 10“® P |9l 

[To]. Long distance contributions overshadow the short distance one with BR{D —)■ 
Xue+e-)^SM - 0(10-^) P[9]. 


3 Exclusive decay modes: D ^ 

Previous studies of the these decays were based on the knowledge of non-leptonic 
weak decays of charm mesons to two light vectors and then vector meson dominance 
was assumed to predict rates for D —I /7 m, or a model of charm mesons as 
heavy mesons accompanied by hidden symmetry approach for the vector mesons as 
done in na. There are also QCD sum rules calculation done by authors of Ref. 
na and more recent one in Ref. na. It was found that the amplitudes fulfil rela¬ 
tions A{D^ —)■ K*^'y) ~ A{D^ —)■ p° 7 ) and A{Df —)■ p''' 7 ) — A{D~^ —)■ p^ 7 ) [I3| 
that led to the predictions for the branching ratios BR{D^ —)■ K*^^) ~ 1.5 x 10“^, 
BR{D^ —)■ p° 7 ) ~ 3.1 X 10“®, BR{Df —)■ ^’*' 7 ) — 2.8 x 10“® and BR(D^ —)■ p^ 7 ) — 
2.7 X 10“®. On the experimental side there only two results for the branching ratios 
BR{D° K*%)exp ^ 3.27(34) x 10"^ and BR{D^ ~ 2.70(35) x 10“^ 

One might update calculation of m including more recent results for the D —V 1 V 2 
helicity amplitudes. However, the relative phases of different contributions are still 
not possible to obtain and only range of the values for the branching ratios can be 
given: BR{D^ —)■ K*^'y) ~ (2.8 —4.9) x 10“^ and BR{D^ —>■ 07) ~ (2.8 —4.1) x 10“^. 
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4 —)► , D —)► Pjji^ jji and D —> P1P2 


The LHCb collaboration improved the bound on the rate BR{D^ —)■ < 

6.2(7.6) X 10“® [15] and for the first time, they determined limits on the branching frac¬ 
tions in several di lepton invariant mass bins in BR{D^ —)• ^~) < 7.3(8.3) x 10“® 

|16j . At the low dilepton invariant mass region 0.25 GeV < < 0.525 GeV the 

LHGb collaboration found upper bound on the rate BR{D^ —)■ pfi~)i.e.b < 2 x 10“®, 
while at high dilepton invariant mass 1.25 GeV < < 2.0 GeV, BRh.e.b{D^ 

< 2.6 X 10“® [I6], at 90% conhdence level. These two results enable to con¬ 
strain size of the Wilson coefficients entering effective Lagrangian ([^. This puts 
then limits on NP contributions in c —)■ ul^l~ in a model independent way. For 
analyses of NP effects in Z1+ —)■ 7 r+/i+/r“ one needs matrix elements of and 

- 1 - 75 ) 0 . We follow here standard parametrisation of these matrix elements de¬ 
scribed in [^[T 8 ]: < 7 r(fc)|M 7 ^(l— 75 )c|Z 1 (p) >= {p+kY+ For the 

f+{qY form factor we use Becirevic -Kaidalov parametrisation [T9| as given in detail 
in [IH]. The tensor current matrix elemet is parametrised as < 'n'{k)\ua^'^c\D{p) >= 
i2fT{qY/'i^D + ^ 7 r[(p + k^YY' “ (P + k-YqY- The HFAG report [^ was used to for 
the relevant parameters present in f+{qY and /t(?^) as given in [18]. Based on the 
effective Lagrangian ([^, the most general expression for the short distance amplitude 
can be written as: 


MsDiD+ip) 7r+(fc)/i+(p+)p (p_)) = 


rrir. 


rriD + rn-,, it 


Cifriq^ 


H— c<^fYq^) 

vr 


u{P-)Pc,pv{p+) + -Ciof+{q ) m ( p _) p „ 7 " 75 ' i ;( p +)} . 

TT 


(3) 


The branching ratio for D —>■ pYp can be written as: 


BR{D^ p+p-) = 


D 


647r3 


\ 


4m^, 2m^, 

1- Y\^C 


m% 


m% 


10 


(4) 


Within SM long distance dynamics can be described be the processes —)■ 

with u and 0 in which then V° decays to p~^p~ pair, presented in details in 

Ref. [IT] for the contribution of D~^ —)■ 7i^p^{u) and updated for the —)■ 7r+0 —)■ 

p'^p~ in Ref. na. The existing experimental upper bound in the non-resonance 
regions indicates that the long distance contribution is fairly suppressed. On Fig. 1 
we present SM contributions to the differential branching fraction for —)• p'^p~ 

as a function of dilepton invariant mass. We also give experimental upper bound for 
the differential branching ratio as found by LHGb [TB] . 

If one considers contributions of NP and if new particle is a new scalar or pseu¬ 
doscalar particle mediating the decay c —)■ then the same particle would con¬ 

tribute to — D oscillations and the physical observable from this process restricts 
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Figure 1: Long distance contributions to differential branching ratio for —)■ 
as a function of dilepton invariant mass (green) and the LHCb bounds 
for the non-resonantt bins (orange). 

the couplings of this operator. The same holds for the flavor changing Z or new 
Z' boson. In the case that NP is generated at the loop level in c —)■ ui^i~ then it 
contributes to — D at the loop level too, as presented in Fig. 2. In addition to dif¬ 
ferential branching ratios at low/high dilepton invariant mass NP detection in these 
decays was also discussed by suggesting new observables. It was found that two an¬ 
gular asymmetries, namely the T-odd di-plane asymmetry and the forward-backward 
dilepton asymmetry offer direct tests of new physics due to tiny SM backgrounds [10] . 

It is important to mention that outside resonance regions of —)■ 7r ~^the 

long distance, as well as SM short distance contributions are more than two orders 
of magnitude smaller than the total branching ratio. Experimental results for the 
differential decay width distribution at the low/high dilepton invariant mass bins 
can be explained by the contributions of the effective Wilson coefficients. We allow 
only one Wilson coefficient at the time to have maximal value. At the same time 
—)■ can give bound on the Cio. Results are presented in Table [Tj It turned 

out that the upper bound on is more restrictive on Ciq Wilson coefficient 

than any of the differential branching ratios for —>■ in dilepton invariant 

mass bins. In Table we use notation Ci = VubV*^,Ci. The differential rate at high 
dilepton invariant mass bin is more restrictive for the effective Wilson coefficients 
than the differential rate at low dilepton invariant mass bin. Many models of NP 
were discussed in literature as supersymmetry with or without R-parity violation, 
new vector-like quarks, leptoquarks imiiH], Little Higgs models [101 . All of these 
models modify some of Wilson coefficients Cj, i = 7,9, 10, however, they are still 
smaller than the bounds given in Table [Tj The detailed analysis of the semileptonic 
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c — ul'^i 


D° -D° 






Figure 2: NP contributions to c —>■ ul^l and 



at tree and loop level. 


four body D —)• hhl^l~ decays was done in the work of Ref. [21]. The dominant 
long distance contributions (bremsstrahlung and hadronic effects) are calculated and 
total branching ratios and the Dalitz plots are presented. Assuming vector meson 
dominance, it was found BR{D^ —)■ ~ 10“®, BR{D^ —)■ ~ 

10-®, BR{D^ K-K+l+R) ~ 10-^ and BR{D° A:+7r-/+/-) ~ 10"®. 


5 Direct CP violation in rare D decays 

In 2011 the LHCb collaboration [5] found rather large CP violation in the difference 
of CP violating asymmetries for D —>■ 7r~^7i~{K~^K~). Using the updated result of 
the LHCb collaboration [22|, HFAG [20] produced the world average CP asymmetry, 
AAcp = AcpiD K+K-) - Acp{D tt+tt-) = (-0.253 ± 0.104)%. Whether CP 
violating asymmetry is present in H —)■ tt+tt" (iC+iC”) or not is still an open issue. 
For charm meson decays, the CP violating asymmetry is dehned as: 

, p(D0 ^ /) 

Acp — -=n-• 

T{D^^ f) + T{D^ ^ f) 

The authors of Ref. |23| investigated CP violating observables in H —)■ Uy —)■ P+P "7 
decays. If CP violation is due to NP effects [23|, then it is most likely a result 
of the chromomagnetic operator Qg contribution m- In the case of rare D de¬ 
cays CP violation results from mixing of Qs into Qi under QCD renormalization. 

{mc)]\ — Im[C^'^(rUc)]I — 0.02 x 10"^. NP contribution can be comparable 
in size with the real part of the SM {m^\ = (0.5 ±0.1) x 10"^. This implies 
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Max. Ci 

{D -)■ 7riJ,n)i,i,b. 

{D 'P^iii)h.i.b. 

D fxfd 

\C7\ 

0.54 

0.46 

- 

\C9\ 

1.33 

0.91 

- 

\Cio\ 

1.32 

0.68 

0.63 

IG 9 I = -|Gio| 

0.91 

0.54 

- 


Table 1: Maximally allowed value of the Wilson coefficients, Ci = VubV^Ci^ calculated 
in the non-resonance regions of D'^ —)■ at the low lepton invariant mass 

^ [0.25,0.525] GeV, denoted by l.i.b. and at the high invariant mass region 
G [1.25,2.0] GeV, denoted by h.i.b., and from the upper bound on the rate of 
Br(D° —)■ < 7.6 X 10“®. The last row gives the maximal value for the case 

IG9I = -|Gio|. 


that if the phase of long distance contribution can be neglected and the relative strong 
phase is maximal, the GP asymmetry can reach the 0(1%) level. The current world 
average of the GP violating asymmetry in AAcp, following the work of [23] leads to 
a GP asymmetry in D ^ K^K~'y of the order 1%. 

We found that in H —)■ there is a possibility to study GP violation ob¬ 
servables [T 8 |. New GP violating effects in rare decays D —)■ are consequence 

of the interference of resonant part of the long distance contribution and the new 
physics affected short distance contribution. The observables, the differential di¬ 
rect GP asymmetry and partial decay width GP asymmetry can be introduced in 
a model independent way. Among all decay modes the simplest one for the exper¬ 
imental searches are D~^ —)■ and Df —)■ Only when third gen¬ 

eration is included there is a possibility to obtain non-vanishing imaginary part: 
Im(A;,/Ad) = —Im(A 5 /A(i). The GP violating parts of the amplitude are suppressed by 
a very small factor \b/\d ~ 10“^ with respect to the GP conserving ones and there¬ 
fore the GP violating effects should be very small. In vicinity of the 0 resonant peak, 
the long distance amplitude for —)■ decay is well approximated by non- 

factorizable contributions of four-quark operators in The width of 0 resonance 
is very narrow (P^/m,^ 4 x 10“^) and well separated from other vector resonances 

in the spectrum of iA —)■ Pi^i~. Relying on vector meson dominance hypothe¬ 
sis the g^-dependence of the decay spectrum close to the resonant peak follows the 
Breit-Wigner shape mm With A'Ijj = the differential direct GP violation 
becomes 


dCP 



|MP + |M|2 




( 6 ) 


The asymmetry can reach acp ^ 


1% (see discussion in mi). 


In addition, a GP 
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Decay mode 

size 

Reference 

D —>■ p{uj)'y 

<3% 

P! 

D K+K--f 

< 1% 

123 

D XJ+l- 

<3% 

[?] 

D~^ — )■ pC pT 

<2% 

m 

D~^ —)■ hhp^p~ 

< 1% 

PH 


Table 2: CP violating asymmetries for charm rare decays, the size and the original 
reference. 


asymmetry of a partial width in the range mi < mu < '^2 can be dehned as: 

r(mi < < m2) - r(mi < < m2) 

< n^) + r(,n. < m„ < m,) ’ 

where P and P denote partial decay widths of D~^ and D~ decays, respectively, to 
7r=*"/i+/i“. Acp can be related to the differential asymmetry acp{y/^) as described 
in [18]. The largest possible asymmetries are of the order few percent. The new 
physics detection in these decay modes was also discussed. It was found that two an¬ 
gular asymmetries, namely the T-odd diplane asymmetry and the forward-backward 
dilepton asymmetry offer direct tests of New Physics due to tiny Standard model 
backgrounds. If supersymmetric and Z'-enhanced scenarios are assumed, and if the 
size of Wilson coefficients Cg and Cio is compatible with the observed CP asymmetry 
in nonleptonic charm decays and flavor constraints, it was found in |2l| that new 
physics effects in —)■ hih 2 l~^l~ might reach the percent level. 

If supersymmetric and Z'-enhanced scenarios are assumed and if the size of Wilson 
coefficients Cg and Cio is compatible with the observed CP asymmetry in nonleptonic 
charm decays and flavor constraints, it was found in |2T] that new physics effects in 
> hih 2 l^l~ might reach the ~ 1% level. 

6 Summary 

Within SM rare charm decays are fully dominated by long distance dynamics. Recent 
results of LHCb experiment on H —)■ and —)• fjCfjT enable to determine 

bounds on effective Wilson coefficients: Cj, i = 7, 9, 10. 

It was found that some signals of new physics might arise in H as well 

as in decays with the leptonic pair in the hnal state D —)• —)■ , 

—)■ In discussion of NP in rare charm decays it is necessary to check 

whether oscillations give additional constraints on the couplings of new 

physics. At the same time, K and B physics might for doublets of up-like quarks 
interacting with NP particle, give very restrictive bounds. 
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The possible presence of CP violation indnced by new physics in charm nonleptonic 
decays has stimulated a number of studies. The three body D —decay decays 
is particularly interesting, since one can focus on the CP asymmetry around the 0 
resonant peak in spectrum of dilepton invariant mass. The appropriate observables, 
the differential direct CP asymmetry and the partial decay width CP asymmetry can 
be introduced in a model independent way. 

Although long distance dynamics overshadows short distance contributions in rare 
charm decays, more precise measurements and improved knowledge of hadronic quan¬ 
tities, might uncover presence of New Physics. 
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